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Abstract. We present a new method to analyze anisotropic flow from the genuine 
I correlation among a large number of particles, focusing on the practical implementation 

of the method. 



1. Introduction 
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O ■ Anisotropic collective flow is defined as a correlation between the azimuthal angle (j) 
O ' of an arbitrary particle and the azimuth $/j of the impact parameter in a non-central 
nucleus-nucleus collision. It is commonly characterized by the Fourier coefficients of the 
single-particle distribution [1] 

O 

5 : Vn = (cos(n(0 - , (1) 

where (. . .) here denotes an average over many particles and events. In particular, 
I elliptic flow V2 is recognized as a sensitive probe of thermalization at RHIC [2]. 
H ' While anisotropic flow is by definition a collective effect that involves many 

particles, it has mostly been analyzed using methods based either on a study of two- 
particle correlations [3-5] or on the cumulants of correlations between a few (in practice, 
up to 8) particles [6]. We recently proposed a new method of analysis [7,8] that 
remedies this limitation, and extracts flow from the correlation between a large number 
of particles instead of only a few. In the following, we introduce the practical recipes for 
implementing the method, referring the reader to Refs. [7,8] for theoretical justifications. 

2. Integrated flow 

For a given centrality bin, the first step in the analysis is to obtain an estimate of the 
flow integrated over some phase-space region (typically corresponding to the acceptance 
of a detector or a set of detectors) . We define integrated flow as the average over events 

Vn = Wj cos(n(0j - , (2) 
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where the sum runs over all particles detected in an event. As is usual in flow analyses, 
Wj is a weight attributed to the j-th particle so as to maximize the integrated flow 
value. Further details regarding the choice of weights can be found in Ref. [6]. With 
unit weights, and neglecting for simplicity multiplicity fluctuations, Vn = Mvn, where f„ 
is an average of the Fourier coefficient ((H) over the phase space covered by the detector. 

To derive an estimate of the integrated flow Vn (for practical purposes, n = 2 at 
RHIC and LHC, n = 1 at AGS and below), one should first compute for each event the 
complex-valued function^ 

M 

g\ir) = Y[[l + irWjCOs{n{(pj - 9))] (3) 
i=i 

for various values of the real positive variable r and of the angle 6 {0 < 6 < n/n; 
in practice, 4 or 5 equally spaced values of 6 are enough as shown in Ref. [8]). The 
(f)j are the measured azimuthal angles of the particles, using a fixed reference in the 
laboratory, and the product runs over all particles. Please note that the method is 
stable against effects, like multiple hits or showering, which bias the results of other 
methods of analysis, so that one should not refrain from using all detected particles, 
combining information from different detectors: increasing the multiplicity results in 
smaller statistical uncertainties on the flow estimates. 

Together with values of g^{ir), one should in each event compute the sums 

M M 

Qx = Yl cos(n0j), Qy = Y^ Wj sm{n(j)j), (4) 
i=i j=i 

as well as their squares Qj^ and Qy'^. 

Next, one should average g^iyr) over events for each value of r and 9: 

J^evts events 

where iVe^ts is the number of events used in the analysis. This is also a good time to 
compute the averages over events {Qx), (Qy) and {Qx"^ + Qy^)- 

For every 9 value, one must then look for the position Tq of the flrst positive 
minimum of the modulus |G^(ir)|. An estimate of the integrated flow Vn is given by 

{oo} = (6) 

where joi — 2.40483 is the flrst zero of the Bessel function Jq. If the detector acceptance 
has reasonable azimuthal symmetry, the estimates do not depend on 9 up to statistical 
fluctuations (see below). § One eventually averages l^{oo} over 9. This yields a new 

X In Refs. [7,8] we used a different generating function. The one defined in Eq. Q, which is introduced 
in Appendix A of Ref. [8], actually yields more accurate results for higher harmonics V2m v^n when 
analyzing differential flow. However, both forms of generating functions are actually equivalent for the 
most part, and the computations of Ref. [8] can easily be adapted to the present function Eq. 
§ For an anisotropic detector, T/f {oo} shows an oscillatory pattern that can be computed in terms of 
the Fourier coefficients of the detector acceptance-efficiency profile. There is an extra proportionality 
factor between y„{cx)} and Vn that can also be calculated and is close to unity in most cases, see Ref. [8] 
for details. 
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estimate y„{oo} with smaller statistical errors. This estimate is then used to compute 
the resolution parameter, which measures the relative strength of flow compared to 
finite-multiplicity fluctuations: x = Ki{oo}/cr [4], with a given by 

(7) 



(8) 



a' = {Q,' + Qy') - {Q^y - {QyY - K{oo}^. 
The relative statistical error on V^{oo} due to the finite number of events is 



2iVe.tsJoi^l(jOl)2 
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where Ji is the spherical Bessel function of order 1. The statistical uncertainty on the 
estimate Vn{oo} is a about a factor of 2 smaller [8]. 



3. Differential flow 



Once integrated flow estimates have been obtained, one can turn to the analysis of 
differential flow, i.e., the flow of particles of a given type in a definite phase-space 
window, which we shall call "protons" for the sake of brevity. A "proton" azimuth 
will be denoted by ifj, and the corresponding differential flow Vp{pT,y) by v'^. Using an 
estimate of integrated flow in harmonic n, as e.g. V^^{oo}, one can derive an estimate 
of Vp in any harmonic p which is a multiple of n, i.e., p = mn with m integer. 

Now, for a given angle 9, with the help of the values of Tq and V^{oo} determined 
following the recipe of Sec. El an estimate of v'^^ is given by 

/ / a,, a. cos{mn{ip — 9)) \ \ 



^{oo} _ Ji{f 



01 J 



Re 



01 > 



o; 
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Wj cos(n 



(0, - 9)) 



(9) 



+ ir^Wj cos{n{(j)j — 9)) / j 



In the denominator, the average is over events and the sum runs over all particles in 
each event. By contrast, the average {■ ■ in the numerator is over protons, and we 
have denoted by the weight associated with a proton. Please note that the sum in 
the denominator need only be computed once per event (it is actually the derivative of 
g^{ir) at the minimum Tq), while the quantity to be averaged in the numerator varies 
from one proton to the other, even for protons within the same event. Finally, Re 
denotes the real part of the (complex-valued) ratio. 

Denoting by A^' the total number of "protons" in the phase-space bin under study, 
the statistical uncertainty on the estimate v'^^^{oo} is 

^ ^expf ^] + (-irexp(-^] Jo(2joi) 



4iV'J^(joi)2 



(10) 



As in the case of integrated flow, averaging the various estimates f'^„{oo} results in a 
new estimate w^„{cxd} with reduced statistical error bars (by a factor ^ 2). Regarding 
the effects of detector anisotropics, they are the same as above too [8] : a ^-dependence of 



f'„„{oo} and a multiplicative factor between the "true" and its estimate f^„{oo}. 
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Finally, let us briefly mention systematic errors inherent to the method. A careful 
study in Ref. [8] allowed us to conclude that the relative error due the interplay of 
nonflow effects and higher harmonics (in particular f2„) is of order 



The relative error on integrated flow estimate V^{oo} is of the same order of magnitude. 
Actually, in the analysis of higher harmonics (m > 1), an extra error term arises 



Unlike the previous term, this is an absolute systematic error, not a relative error on 
the flow. With any other method of flow analysis, the systematic error will always be 
larger (or at least equal). 

4. Discussion 

The method of analysis we presented above is simple to implement: compute values of 
the generating function (0) and find the first minimum of its modulus, then Eq. (jHI) gives 
you the integrated flow. Knowing the position of the minimum, you can then compute 
in a second pass through the data the quantities in the right-hand side of Eq. (jH)) and 
deduce differential flow. Nothing more is required. 

In addition, the method is conceptually rich. The minimum of |G^(ir)| is in fact 
compatible with a zero of G^{z), where z is a complex variable. The mere existence 
of such a zero close to the origin and scaling with the inverse of the system size {Vn in 
Eq. is roughly proportional to the multiplicity M) signals the presence of collective 
effects. This is analogous to Lee- Yang theory of phase transitions [9] in which the zeroes 
of the partition function come closer to the origin with increasing system size if, and 
only if, there is a phase transition. Since anisotropic flow is a collective effect, Lee- Yang 
zeroes definitely are the most natural method to analyze flow. 
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